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ABSTRACT 

With the use of a background Milky-Way-like potential model, we 
performed stellar orbital and magnetohydro dynamic (MHD) sim¬ 
ulations. As a first experiment, we studied the gaseous response 
to a bisymmetric spiral arm potential: the widely employed co¬ 
sine potential model and a self-gravitating tridimensional density 
distribution based model called PERL AS. Important differences 
are noticeable in these simulations, while the simplified cosine po¬ 
tential produces two spiral arms for all cases, the more realis¬ 
tic density based model produces a response of four spiral arms 
on the gaseous disk, except for weak arms -i.e. close to the lin¬ 
ear regime- where a two-armed structure is formed. In order to 
compare the stellar and gas response to the spiral arms, we have 
also included a detailed periodic orbit study and explored different 
structural parameters within observational uncertainties. The four 
armed response has been explained as the result of ultra harmonic 
resonances, or as shocks with the massive bisymmetric spiral struc¬ 
ture, among other. From the results of this work, and comparing 
the stellar and gaseous responses, we tracked down an alternative 
explanation to the formation of branches, based only on the orbital 
response to a self-gravitating spiral arms model. The presence of 
features such as branches, might be an indication of transiency of 
the arms. 

Key words: MHD - Galaxy: disc - Galaxy: strcuture - Galaxy: 
kinematics and dynamics - galaxies: spiral - galaxies: structure 


Spiral arms are one of the the most striking, beau¬ 
tiful and scientifically challenging structures of 
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1 INTRODUCTION 


disc galaxies. They have fascinated and intrigued 
astronomers for centuries. As a rough initial ap¬ 
proximation to treat this complex problem, the 
spiral arms had to be assumed as almost mass¬ 
less and/or with extremely modest pitch angles, 
in order to obtain a solution as a linear perturba¬ 
tion to the axisymmetric background potential. 
The proposed solution to the permanence of spi- 
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ral arms in galaxies wa s based on the sp iral den¬ 
sity wave linear theory (|Lin Sv Shu II 1964 inspired 
in the 60’s work of B. Lindblad and P.O. Lind- 
blad). The analytical solution at first order of 
the theory, known as the tight-winding approx¬ 
imation (TWA), that represents a weak poten¬ 
tial model, extremely idealized as a smooth, even 
negligible perturbation to the background poten¬ 
tial, was modeled as a simplified cosine function 
to represent the spiral arms gravitational poten¬ 
tial. Probably followed by this initial attempt of 
solving the nature of spiral arms, a lot of work 
has tackled repeatedly the problem of modeling 
them, as a simplihed periodic function that disre¬ 
gards the importance of their dynamical effects. 
However, spiral arms have proven to be an influ¬ 
ential feature on galactic modern models (such 
as PERLAS and N-body models, based on three 
dimensional density structures), far beyond of a 
simple smooth perturbation. 


If we compare the spiral arms with the galac¬ 
tic bar, the latter has a mass between 10 and 
20% of the disc mass (iM atsumoto et al.l Il982l: 
Dwek et al.lll995l: IZhaol[l996l : IWeiner Selwood I 


1999|), unlike the spiral arms that have a mass 
smaller in general th an 5% of the disc mass 
(|Pichardo et al. ir2003l l . Consequently, the largest 
influence on the disc of strongly barred galaxies 
is, of course, due to the central bar. Therefore 
one would be tempted to oversimplify the sce¬ 
nario, by assuming that spiral arms are not influ¬ 
ential at all on the dynamics of the disc. However, 
there are determinations that suggest that the 
majority of spiral galaxies, or at least those with 


beini 

y linear (lAntoia et al. 1 1201 ll: lAntoia et al.l 

20ld 

: Antoia et al. 

2009 

: IVofflis et al. l2006l: 

Roca-Fabrega et al. 

201.1 

2014: IPichardo et al. 

2003 

: iKawata et al. 

|2014 ISellwood & Carlbere 

2^ 

: ISellwoodll2011 



Based on all those rigorous studies, spiral 
arms seem to merit an extra effort to model them 
beyond a simple perturbing term. This fact is 
even more significant in the case of the inter¬ 
stellar gas, which is considerably more sensitive 
to the details of the potential, and in general, 
responds strongly even if the mass of the spiral 
arms is only a small fraction of the axisymmetric 
background (iKim et al.ll20l3 : 1 Gomez et al .11 201 
[SM et al ][i92T 


For almost a century, formal studies 
of spiral arms have been carried on. Yet, 
the spiral arms morphology, origin and na¬ 


ture are still poorly understood. With re¬ 
spect to their morphology, there is plenti¬ 
ful theoretical and observational literature 
on structures related to spiral arms, such 


and their plausible explanations (IWeaver 

1970 

: Lvndd 1970l: Shu et al.l 1973: 

Sofue 

1976 

; lElmegreen Il980l; IScoville et al.l 

200l|; 

Lee Shu l2012l; Kim Ostriker 1 

20061: 

Chakrabarti et al. 12003 

; iKim Ostriker 

2002 

tiDobbs et al.ll201ll; 

Shettv Ostriker 

200d 

1; ICorder et al.l 20081; iBalbusI Il988ll. 


Regarding spurs and feathering, interest¬ 
ing scenarios to explain these features 
include plain hydrodynamic simulations, 
where the Kelvin-Helmholtz instability 
creates spurs pr oduced by shocks w ith 
the spiral arms (jWada &: Kodal l2004l ; an 
extension to self-gravitating three dimen- 
sio nal models of this wo rk is presented 
in iKim Ostrikerl l20Q6l T Other scenar¬ 
ios include the use of global MHD simu¬ 
lations, adding gas self-gravity and mag¬ 
netic fields that produce differential com¬ 
pression of gas flowing through the arms 
resulting in the formation of sheared struc¬ 
tures in the interarm r egions that resem¬ 
ble s p urs and feathers (IShettv Ostrikerl 
l2Q06l ; iKim Ostrikerl l2002h . or through 
grav itational or magneto-Jeans instabili¬ 
ties (|Balbuslll988l T 

On the other hand, features such as 
branches are significantly longer than 
spurs that emerge almost parallel to the 
main spiral arms but with different pitch 
angles (usually smaller), and are generally 
associated with resonances. It is worth 
to note that, in the literature, the names 
spurs, feathers and branches are assigned 
to slightly different structures. In this 
work we will adopt the definiti on given 


by Feitzinger Schwer d tfegerl (|l982l l 


and Ichakrabarti et al. I ([200^ 1 . where the 
“spurs” and “feathers” are small struc¬ 
tures that arise almost perpendicularly 
from the main spiral arm, and “branches” 
are narrow structures that arise almost 
parallel to the main spiral arm, with 
smaller pitch angles and extension longer 
than spurs (examples of branches can be 
seen in galaxies like NGC 309, NGC 1637, 
NGC 2997, NGC 6946, among others). 
Regarding such structures, theory and 
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gas simulations have already shown a 
bifurcation of spiral arms produced by 
ultraharmonic resonances induced by 
the main s piral arm on th e background 


the mam s piral arm on th e background 
gas flo w (fshu et al] 19731; Patsis et alJ 


1994 [T 997 I: IChakrabarti et al. I l2003l: 


Artvmowicz Lubow 19921 ^. In this case 


the bifurcation of the arms is expected 
due to the topology of the stable peri¬ 
odic orbits at the 4/1 resonance. Other 
interpretations also include the possibility 
that branches are the response of the 
gas produced by strong shocks against 
the main massive spiral arms, i.e. the 


stars (iFuiimotol 19681; RobertsI 

1969; 

Shu et aLlll973l; iBarbanis Woltierl 

1967; 

Roberts, Huntlev van Albada 1 

1979; 

Martos et al.l l2004l; iKim Kim 1 

2014; 

Chakrabarti et al. 1 l2003l; lYanez et al. 

20081). 


Specifically, regarding the Milky Way 
Galaxy, its general structure has been exten¬ 
sively studied and the latest determinations 
seem to agree with two grand design symmet¬ 
ric spiral arms, seen on IR and optical, and 
some other weaker arms, detectable mainly 

in _ op tical observatio ns _ ( Drimmel &: Sper^ 

200 ll: iB eniami n et alJ [2005: IPrimmel l2000l : 


Drimmel Soergell l200ll : IValled 120131 . l2002l b 
The weaker arms (branches) have been ex¬ 
plained through hydrodynamic simulations as 
the response to the two-armed stel l ar pa ttern. 
For e xamp le, lEnglmaier GerhardI (|l999l ) and 
IFuxI (|l999l ) found that the gas response to a 
barred potential can produce a four-armed spiral 
structure. Likewise, from hydrodynamic and 
magnetohydrodynamic (MHD) simulations, the 
gas responds to an imposed two- armed spiral 
poten t ial with four spira l arms ([Gomez et al.l 
1201.1 iMartos et af] 120041 : IShu et alJ Il97.1 . In 


this context, the gas component shocks at the 
position of the spiral arms, prod ucing a density 
enhancement a n d star format ion imotoll 19681: 

RobertsI Il969l : iMoore et al.l I 2 OI 2 I : Seo Kinj 


2 OI 4 I I that, depending on the relative velocity 
between the gas and the spiral arm, dwell 
upstream or downstream or in the spiral arm. 
In particular, iRobertsI (| 19691 ) showed that the 
nonlinear response of the gas to a stationary 
stellar spiral arm potential may produce two 
shocks, which then would be associated to a 
doubling of the spiral arm number as seen in the 


gas component. IShu et aTI (|l973l ) demonstrated 
that shocks in galaxies arise necessarily if the 
spiral arm strength exceeds a certain critical 
value. Additionally, they found a range of val¬ 
ues for the wave frequency that generates an 
ultraharmonic resonance which can provoke a 
secondary compression of the interstellar gas. 
This effect has been related to the origin of th e 
Garina arm in the Milky Way (|Shu et al.lll973l ). 
for example. However, since branches have been 
traditionally explained as shocks induced by the 
spiral arms, it is important to mention other 
observational work that shows there is little 
difference (or non at all) in the star formation 
efficiency (i.e. no shocks) betwe en the spiral arms 
regions and the r e st of the disk (iFoyle et al.ll201^: 
Eden et al I 2 OI 3 I : IEovIc et ffi] 1 20 111 : I Dobbs et all 
201 ih . 


Einally, regarding their nature, a prob¬ 
lem of great interest is whether spiral arms 
are a transient or a long-lived feature. On 
this matter, recent numerical simulations show 
that spiral arms are transient a n d recurrent 


201 ll: Roskar et al.l h 

2 OI 2 I: iD’Onghia et al.l 

2013: 

Perez-Villeeas et al.l 

201 I I 2 OI 2 I: ISellwoodI 

2011: 

Sellwood & Garlbers 

20 id: Kawata et al.l 

2014; 

Eovle et al.l I 2 OII 

). However, in a stud 

1 

Scarano & Leoind (|2013l) on a samole of 27 { 


ies, the authors concluded that the break found 
in the radial metallicity distribution near to coro¬ 
tation resonance (GR), implies that the spiral 
structure is a rather long lived feature. Other 
authors, employing different techniques, seem 
to find o bservational proofs of long-lasting spi¬ 
ral a r ms (iMartm ez-Garcfa Sz G onzalez-Lopezliral 
I 2 OI 3 I : IPonner Thomassonlll994l : IZhandll998l ) 


In this work we explore the gas response to 
a bisymmetric spiral arm potential. Eor this pur¬ 
pose we employed two different models, a cosine 
potential (generally used in literature) and the 
three-dimensional, density distribution based po¬ 
tential PERLAS, applied to the particular case 
of a Milky Way-like galaxy. With this study we 
find an alternative explanation to the formation 
of branches in disk galaxies, and their relation to 
transient or long-lived spiral arms. 


This paper is organized as follows. In Sec¬ 
tion [2l the galactic potential used to compute the 
stellar orbits and the MHD initial simulation set¬ 
up are described. In Section O we present first 
a comparison between the models of the spiral 
arms: PERLAS and the cosine potential; second. 
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4 Perez-Villegas, Gomez & Pichardo 


we present the gas response to spiral arms PER- 
LAS model changing the structural parameters 
such as the pitch angle and the mass of the spiral 
arms, and their connection with the presence of 
galactic branches as a signature of transient spi¬ 
ral arms. Finally, we present a discussion and our 
conclusions in Section |4l 


shapes the starting point of the spiral arms, so 
that N ^ 0 represents a 180° tip transition from 
the bar to the spiral arms, while N ^ oo means a 
90° tip transition from the bar to the spiral arms. 
In the present wor k, we set it equal to 100 (for 
further details, see lPichardo et al7ll2003h . 


2 METHODOLOGY AND 

NUMERICAL IMPLEMENTATION 

Motivated by the fact that the gas in galaxies is 
dynamically colder than the stellar disk, in addi¬ 
tion to being collisional, we can expect it to be 
much more sensitive to details of the given po¬ 
tential than stars. Therefore, any differences be¬ 
tween the cosine and PERLAS potentials might 
be magnified in gas with respect to the stellar re¬ 
sponse. And so, for purposes of comparison, we 
produce the whole study employing both poten¬ 
tials. In this section we introduce brieffy the co¬ 
sine and PERLAS potential models and the MHD 
setup applied to a Milky Way like Galaxy. 

In all cases, the spiral arm models are super- 
impos ed to the axisymme t ric ba ckground poten¬ 
tial of I Allen Santillan I (|l99lh . which includes 
a Miyamoto-Nagai bulge and disk, and a super- 
massive spherical halo. Table [T] presents the basic 
parameters of the axisymmetric background po¬ 
tential. 


2.1 Spiral arm models 

2.1.1 Cosine Potential 

As mentioned before, a large majority of the in¬ 
vestigations on the gas associated to spiral arms 
model them as a linear perturbation of the ax¬ 
isymmetric background, represented by 

^sp{R, 0) = f{R) cos [20 + g{R)] , (1) 

where R, 0 are cylindrical coordinates, f{R) is 
the amplitude function of th e pertu rbation, given 
by IContoDOulos Grosboll (|l986h as f{R) = 
—ARe^^^, where A is the amplitude and is the 
inverse of the scale-length. Finally, g{R) describes 
th e geometry of the spiral pattern ( l ocus), given 
by iRoberts. Huntley van Albada 1 (I 1979 I 1 , as 

where ip is the pitch angle, Rs is the start posi¬ 
tion for the spiral arms, and A is a constant that 


2.1.2 PERLAS model 


In contrast, the PER LAS spiral arms model 
(| Pichardo et al. I l2003l l is a bisymmetric, three- 
dimensional, self-gravitating stationary model 
based on an adjustable mass distribution, rather 
than a local arm approximation, as the co¬ 
sine potential. Several studies have shown that 
there are differences in the stellar orbital dy¬ 
namics when the spiral arms are modeled with 
the cosine potential or with the PERLAS mode l 

(|Pichardo et al. |[2003l : lAntoia et al. II2OO9I . l20Il| j . 

Furthermore, chaotic orbital studies have demon¬ 
strated that a more detailed spiral arm potential 
(e.g.,the PERLAS model) induces an important 
fraction o f chaos, enough to destroy the spiral 
structure (|Perez-Villegas et al.ll20l3 . l20I2l j , when 
all chaotic behavior was originally attributed 
to effects produced by the bar, such as over¬ 
lapping of resonances ([Contopoulos Harsoulal 
I2OI2I : IContoDOulos III967 I l 

PERLAS is constructed by a superposition of 
individual oblate inhomogeneous spheroids along 
a given locus (the same than in the cosine poten¬ 
tial, equation!^ and superimposed to the axisym¬ 
metric background. In this model, the spiral arms 
have a well defined mass, unlike the cosine poten¬ 
tial, where the spiral arms are treated as a peri¬ 
odic function of the potential not straightforward 
translatable to mass. The mass assigned to build 
the spiral arms in the PERLAS model is sub¬ 
tracted from the disk mass. Thus, the inclusion of 
PERLAS to the Galactic model does not modify 
the total mass of the original axisymmetric back¬ 
ground. In Table [ 1 ] we present the parameters of 
spiral arms that we used in our simulations. For 
further details a bout the PERLAS mod el, we re¬ 
fer the reader to I Pichardo et al. I (|2003l j . 


2.2 Force fitting 

In order to compare PERLAS with the cosine po¬ 
tential, we need to fit the spiral arm strength to 
make them as similar as possible. We achieve this 
by adjusting the amplitude of the cosine poten¬ 
tial (the factor A in eq. [T]) so that the resulting 
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Table 1. Parameters of the Galactic Model 


Parameter 

Value 

Reference 

Axisymmetric Components 

Ro 

8.5 kpc 

I 

00 

220 km s“^ 

I 

Bulge mass 

I.4I X lO^o M© 

I 

Disk mass 

8.56 X IQio M© 

I 

Halo mass 

8.002 X 10^1 M© 

I 

Spiral Arms 

Locus 

Logarithmic 

2, 3 

Arms number 

2 

4 

Pitch angle(z) 

15.5 ±3.5° 

4 

Alarms / FI disk 

0.03 ±0.02 

3 

Mass 

2.7-5.4 X 10^ M© 

3 

Inner limit 

3.3 kpc 

ILR position based 

Outer limit 

12 kpc 

GR position based 

Scale-length 

2.5 kpc 

Disk based 

Pattern speed (Dp) 

-20km s“^ kpc“^ 

5 


Referenc es —(1) |AlleiL.fc _SaiitillaijJ [l99l|; ,i2l | Grosb|z)l^ fc Patsis 


^98; IPichardo et alTl20^ : fd forimmel II2OO0I : iMarto^^d 


2004. 


radial and azimuthal forces approximate the ones 
obtained using the PERLAS model. 

Figures [Hand [ 2 ] show the resulting fit for the 
case of a Milky Way-like galaxy, where the pitch 
angle is 15.5° and spiral arm mass (in PERLAS) 
is 3% of the disk mass. In this case, the cosine 
amplitude A is 650 km^ sec“^ kpc“^. In Figured] 
each panel represents different radial lines, start¬ 
ing in 0° up to 75°. In the Figure [2| each panel 
shows the computed azimuthal force at different 
radii, from 4 to lOkpc. Continuous lines repre¬ 
sent PERLAS, dotted lines represent our fit to 
the cosine potential. 


2.3 Spiral arm strength 


The spiral strength is related to the pitch an¬ 
gle and the mass of the arm. In the PERLAS 
model, the spiral arm mass is a small fraction of 
the disk mass. To quantify the strength of the 
spir al arms, we calculated the Qt(R) parame¬ 
ter ([Sanders Tubbs |[l980l : ICombes Sanders I 
[l98l|), frequently used to quantify the strength 
of bars and spiral arms (|Buta Block I 1200 ll : 
iLaurikainen Sal^ l2002l i . The value of Qt{R) 
is given by 


Qt(R) 


%R) 


liFniRW 


(3) 


where = | {d^{R,0)/d0) /R|max represents 

the maximum amplitude of the tangential force 
at radius R, and {Fr{R)), is the average axisym- 
metric radial force. Figure [3] shows Qt{R) for the 
PERLAS (solid line) and cosine (dotted line) po¬ 
tential models, for the case where the spiral arms 
mass of the PERLAS model has a 3% of the disk 
mass. 


As opposed to Figures [T] and [2l where the 
force amplitude is almost the same. Figure [3] 
shows that there is a difference between both 
potentials. While the maximum value of Qt for 
PERLAS is ~ 0.096, for cosine is only ~ 0.035. 
Therefore, if the arm strength is measured using 
equation © instead of the arm force, it is nec¬ 
essary to increase the cosine model amplitude. 
Figure m shows the Qt parameter corresponding 
to an increased cosine amplitude, in this case, 
A = 2000 km^ sec“^ kpc“^. In this case, since Qt 
for the cosine model is larger than PERLAS along 
all the radial range, this case should be considered 
as an example to test the gaseous disk response 
to an extreme cosine model, in order to verify 
if at this forcing the cosine can reproduce what 
PERLAS does. In this way we are bracketing the 
cosine potential within the values of the force for 
the PERLAS model. 
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6 Perez-Villegas, Gomez & Pichardo 



R (kpc) 

Figure 1. Radial force of the spiral arms along of the different radial lines, starting in 0° up to 165°, each 15°. 
The radial force is scaled to the axisymmetric background radial force at the solar position. Continuous lines 
represent PERLAS, dotted lines represent the best fit with equation ^ 


2.4 The MHD set-up 

The initial set-up of the simulations consisted 
on a gaseous disk, with density profile given by 
n{r) = noexp[—(r — ro)/rd], where no = 1.1 cm“^, 
ro = 8 kpc and Pd = 15 kpc. The gas follows 
an isothermal equation of state with tempera¬ 
ture T = 8000 °K. Additionally, the gas is per¬ 
meated by a magnetic field, initially in the az¬ 
imuthal direction, with an intensity given by 
B{r) = Boexp[—(r — ro)/rB], where Bq = 5/xG 


and tb = 25 kpc. The disk is set up in rota¬ 
tional equilibrium between the centrifugal force, 
the thermal and magnetic pressures, magnetic 
tensi on and the background a xisymmetric poten¬ 
tial ([Alien Santilld^ Il99lh . This equilibrium 
is perturbed by the spiral arm potentials un¬ 
der study, both rotating with a pattern speed 
Qp = 20 km s“^ kpc“^. 

We employed the ZEUS code 
([stone Sv Normanl Il992al li3) to solve the MHD 
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Figure 2. Azimuthal force at the different radius, from 4 to lOkpc, each 2kpc. The radial force is scaled to 
the axisymmetric background radial force at the solar position. Continuous lines represent PERL AS, dotted lines 
represent our fit with equation ^ 


equations, which is a finite difference, time 
explicit, operator split, Eulerian code for ideal 
MHD. We used a 2D grid in cylindrical geometry, 
with R G [1.5,22]kpc and a full circle in the 
azimuthal coordinate, 0, using 750 x 1500 grid 
points. Both boundary conditions in the radial 
direction were outflowing. All calculations are 
performed in the reference frame of the spiral 
arms. No self-gravity of the gas was considered. 


3 RESULTS 

We present in this section a comparison between 
the gas response to the simplified cosine poten¬ 
tial for the spiral arms, and the density distribu¬ 
tion based model PERLAS. In the limit for weak 
and/or small pitch angles (approximately linear 
regime), both models behave very similar as ex¬ 
pected, however, for stronger arms (more massive 
or larger pitch angles), from this comparison we 
find severe differences in the gas behavior and for- 
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Figure 3. Qt{R) parameter for the spiral arms of a 
Milky Way-like galaxy for the PERLAS (solid line) 
and cosine (dotted line) potential models. 



R (kpc) 


Figure 4. Similar to fig. [S] with the amplitude of the 
cosine model increased. 


mat ion of spiral arms and branch-like structures. 
The deviation of the response between both mod¬ 
els is induced by the basic differences of these po¬ 
tentials. We then present an interpretation on the 
presence of galactic branches prior to the ultra¬ 


harmonic 4:1 resonance and connect them to a 
possible signature of transient spiral arms. 


3.1 Gas Response Comparison: PERLAS 
vs. cosine Potential Models 

In order to compare both potentials, we fit the 
cosine potential with PERLAS to make it as sim¬ 
ilar as possible in the force amplitude f ^2.2l) . In 
ah cases, the axisymmetric gaseous disk, initially 
in rotational equilibrium, was perturbed by the 
spiral potential either PERLAS or the cosine. 

As mentioned already, the cosine poten¬ 
tial represents a simple solution from the den¬ 
sity wave linear theory, self-consistent for tightly 
wound spiral arms (TWA, i.e. where the pertur¬ 
bation is very small, which means small pitch an¬ 
gles or with very reduced masses). The Milky Way 
Galaxy and the most of spiral galaxies are ac¬ 
tually far from this regime. Thus, it should not 
be su rprising that the gaseo us disk (I Gomez et al" 
20^) and the stellar orbits (|Perez-Vihegas et al. 
20121 ) show a different structure when subjected 
to a self-gravitating, more realistic model instead 
of a local approximation. With this in mind, the 
gaseous disk response to both models should be 
similar if we focus on a region of the parameter 
space where both potentials are valid, i.e., if we 
set PERLAS and the cosine potential so that it 
is approximately in the linear regime, with a very 
small pitch angle and spiral arm mass. In Figure 
[5] we show the linear regime for both potentials. 
The density distribution with the cosine potential 
is presented in the left panel and PERLAS model 
in the right panel, the pitch angle is 6°, the spiral 
arms mass is 1% of the disk mass, A = 100 km^ 
s“^ kpc“^. Indeed, the gas response to the po¬ 
tentials is similar as expected, forming two spiral 
arms in both cases. 

For the specific case of a Milky Way-like 
galaxy, we constructed a model that reproduces 
some of the observational parameters for the 
background and spiral arms potential to compare 
them with the cosine potential for the spiral arms 
(force fitted with PERLAS). The spiral arms pat¬ 
tern angular speed is 20 km s~^ kpc“^, on a log¬ 
arithmic locus with a pitch angle of 15.5°. The 
mass of spiral arm for PERLAS model is 3% of 
the disk mass, which is equivalent to a cosine am¬ 
plitude of A = 650km^sec“^ kpc“^, when the 
non-axisymmetric force is employed for the fit¬ 
ting (see ^2.21) . and A = 2000km^sec“^ kpc“^, 
when the arm strength (i.e. the Qt parameter) is 
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X [kpc] X [kpc] 


Figure 5. Density distribution for the simulation with the cosine potential (left panel) and with the PERL AS 
model in the linear regime {right panel) after 1.2 Gyr. In both cases, the pitch angle is 6°. The amplitude for the 
cosine potential is A = 100 km^ sec“^ kpc“^ and the spiral arm mass for the PERLAS model is 1% of the stellar 
disk mass. 


considered for the fitting instead (see ^2.31) . We 
follow the evolution of the system for 5 Gyr. 

Figure [6] shows the resulting density distri¬ 
bution of the gaseous disk when is perturbed 
by the cosine potential with a small amplitude 
{A = 650km^sec“^ kpc“^). After the simulation 
starts, the gas very rapidly settles into a spiral 
pattern, forming two spiral arms at 30 Myr (top 
left panel). 300 Myr into the simulation (top right 
panel), the two spiral arms are better defined, and 
the gas is forming a secondary structure. As simu¬ 
lation progresses (1.5 and 3 Gyr, bottom panels), 
the gas continues responding to the spiral arms 
and even more substructure forms, but the large 
scale density response consists of two arms only. 
The pitch angle of the formed arms is ~ 15.5°, 
equal to the imposed potential at the beginning of 
the simulation. In the last panel, a gas instability 
at the c orotation radius is seen, as reported pr e- 
viously (jGomez et al.|[201^ : iMartos et al.ll2004l i. 

In contrast. Figure 0 shows the density distri¬ 
bution of the gaseous disk when is perturbed by 
the PERLAS spiral arm potential. The gas very 
rapidly settles into a spiral pattern, as in the co¬ 
sine potential case, forming two spiral arms at 30 
Myr. At 300 Myr and later, the gas forms four 
spiral arms instead of the two arms in the cosine 
potential simulation. Even though the simulation 
develops four spiral arms, these are associated in 
two pairs, each composed of a strong arm and 
a weak arm. The strong arms have a pitch angle 
of ~ 15°, and the weak arms have a pitch angle 


of ~ 7°. This doubling of the spiral arms in the 
gas response has been seen in othe r MHD sim¬ 
ulations using the PERL AS model (jGomez et al.! 
l2niA iMartos et, al.ll2nn4l: see also MM- 

In 112.31 we noted that the spiral arm 
strength, as measured by eq. was very dif¬ 
ferent for both potentials, even if the force ampli¬ 
tudes were almost the same. Therefore, we de¬ 
cided to perform a second experiment with a 
larger amplitude A = 2000km^sec“^ kpc“^ for 
the cosine potential, that would be equivalent to 
a factor of ~ 3 times the force of PERLAS. This 
much larger amplitude for the cosine potential is 
stronger than PERLAS at any radii, as shown in 
Eigure [H The purpose of this overestimated ex¬ 
periment was to test if the four spiral arms formed 
with PERLAS were due to the strength of spiral 
arms only. Eigure [8] shows the gas density dis¬ 
tribution resulting from perturbing the disk with 
the larger force amplitude of the cosine poten¬ 
tial. It is readily seen that the effect of the spi¬ 
ral potential on the gas disk is larger and gener¬ 
ates more substructure, but the gas still responds 
forming two spiral arms, unlike the density based 
potential PERLAS that forms four arms for larger 
pitch angles, i.e. for stronger arms. Gonsequently, 
this difference does not come from the force am¬ 
plitude, but it seems rather originated from the 
self-gravitating nature and specific details of the 
potential that a local approximation for the forces 
given by the cosine potential is unable to repro¬ 
duce. We further explain this in the next section. 
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Figure 6. Density distribution using a cosine potential to model spiral arms. Times shown in this plot correspond 
to 30Myr after the simulation starts (top left panel)^ 300 Myr (top right panel), 1.5 Gyr (bottom left), and 3 Gyr 
(bottom right). 


3.2 Branch Formation and their 

Relation with the Transient Nature 
of Spiral Arms 

In ^3.11 we showed that the intrinsic differences 
in the nature of the cosine and PERLAS poten¬ 
tials induce a distinct gaseous arm structure as 
compared to a disk perturbed by a cosine arm 
potential. While employing the spiral arm cosine 
potential, a bisymmetric gaseous structure seems 
to be an invariable outcome, for the PERLAS po¬ 
tential a two or four armed structure could be 
obtained. This result points to the necessity of 
understanding how the gas responds to the PER¬ 
LAS model when we vary the structural parame¬ 
ters of the spiral arms, such as its mass and pitch 
angle, considering the uncertainties in the deter¬ 
mination of these parameters for the Milky Way’s 


spiral arms. Therefore, we developed a set of sim¬ 
ulations varying spiral arm pitch angles (z) and 
spiral arms masses. We present here some of them 
ranging from 12° to 19°, and the mass of the spi¬ 
ral arms (Mgp), from 1.75% to 5% of the total 
disk mass (Md)- 

Eigure[9]shows a mosaic of simulations. The 
panels in the left column are the stable periodic 
stellar orbits and the panels in the right column 
are the gas density distribution after 1.2 Gyr. 

The maximum density response (where 
the orbits crowd producing a density en¬ 
hancement) quantifies the orbital support 
to a given spiral perturbation through 
periodic orbits. We computed the stellar 
periodic orbits in order to explore the 
orbital support to the imposed PERLAS 
spiral arm potential (open squares in 
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Figure 7. Density distribution using the PERLAS spiral arm model. Panels correspond to the same evolution 
times as in fig. [6] 
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the left column-panels of Figure [9|). To 
estimate the dens ity response, we employ 
the m ethod of IContopoulos &: Grosb0ll 
(|l986l ). This method assumes that stars 
in circular orbits in an axisymmetric 
potential, with the same sense of rotation 
of the spiral arms, are trapped around the 
corresponding periodic orbit in the pres¬ 
ence of the spiral arms. For this purpose, 
we calculated a set of central periodic 
orbits (between 50 and 60) and found 
the density response along them using 
the conservation of the mass flux between 
any two successive periodic orbits. With 
this information we seek the position of 
the maximum density response along each 
periodic orbit (fllled squares in the left 
column-panels of flg. [9|). These positions 
are compared with the imposed spiral 


arms (PERLAS model). The method im¬ 
plicitly considers a small dispersion (with 
respect to the central periodic orbit) since 
it studies a region where the flux is con¬ 
served. On the other hand, this dispersion 
is based on parameters for the galaxies 
where dynamics is quite ordered, orbits 
follow their periodic orbit closely, in such a 
way that we consider this study a good ap¬ 
proximation. This method to estimate the 
density res ponse has been widely used in 


literature dContoooulos 8^ Grosbpll 

19881: 

Amaral Leoii 

tiel I 1997 

; Yano et al.l 

20031: 

Pichardo et al. 

bod; 

Voelis et al. 1 

2000: 

Tsoutsis et al.l 

20081; Perez-Villeeas et al.l 

20121. 2013I: 

ljunaueira et al.l l2013lL 


We refer the reader t o th e work of 
IContopoulos &: Grosb^ (|l986l ) for more 
details. With this in mind, a model in 
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Figure 8. Gas density distribution for the simulation with the cosine potential when its force amplitude is 
increased by a factor of ~ 3 (up to the value provided by PERLAS model). The times shown correspond to 30 
Myr (top left panel), 300 Myr (top right panel), 900 Myr (bottom left panel), and 1.5 Gyr (bottom right panel). 


Figure [9] where the open and filled squares 
coincide would represent a stable, approximately 
orbit ally self-consistent system, while a lack 
of coincidence would mea n that the spiral is 
unlikely to be long-lasting (IPerez-Villegas et al.l 

For Mgrp — O.OITSMd and i — 12° (upper left 
panel of Figure [9|, the stellar density response fol¬ 
lows approximately the imposed spiral potential 
prior to the 4/1 resonance. After that, the stellar 
response forms a slightly smaller pitch angle than 
the imposed spiral. In the MHD simulation, the 
gaseous disk responds to the two-arm potential 
with the now familiar four spiral armed structure, 
where the pitch angle of the stronger pair of arms 
corresponds to that of the imposed pitch angle 
potential, while the other pair of gaseous arms 
has a systematically smaller pitch angle, corre¬ 


sponding closely to the regions of periodic orbits 
crowding. Now, with the same spiral arm mass 
but a pitch angle of z = 19° (second row of hg.[9|), 
the stellar density response (i.e. periodic orbits) 
forms again a smaller pitch angle compared with 
the imposed, while in the gas, the four spiral arms 
seem stronger and the difference in the pitch angle 
between the gaseous and imposed arms is larger 
than the previous case. For Mgp — 0.05Md and 
i = 12° (third row), the stellar density response 
closely follows the imposed spiral arm potential 
prior to the 4/1 resonance. After that, the stel¬ 
lar response forms a slightly smaller pitch angle 
than the imposed spiral, while the gas responds 
with four spiral arms, but the second pair is very 
weak but with a significantly smaller pitch angle 
than the imposed spiral arm potential. Finally, 
with the same spiral arm mass (5% of the disk 
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mass) but a pitch angle of z = 19° (fourth row), 
we found no periodic orbits beyond 4/1 resonance 
and the stellar density response forms a pitch an¬ 
gle smaller than 19°. The gaseous disk responds 
with well defined four spiral arms that extend 
up to the corotation radius. Notice that in this 
last simulation, there is not much that can be 
said about the (stellar or gaseous) orbital sup¬ 
port since periodic orbits tend to disappear due 
to the strong forcing of the imposed spiral arms, 
meaning that spiral arms would rather be tran¬ 
sient by construction in this case and the MHD 
gaseous disk behavior is difficult to predict from 
periodic orbit computations. However, such as the 
stellar arms constructed in this case, the gaseous 
spiral would be transient in likely even shorter 
timescales than in the case where periodic orbits 
exist but settle down systematically in smaller 
pitch angles than the original imposed spiral arms 
in the region where periodic orbits do exist. No¬ 
tice that, in general, the arms that should even¬ 
tually disappear in this scenario, are the stronger 
st ellar imposed spiral arms ( see second row of fig. 
[9l IPerez- Villegas et al. I l20ia) . 

Therefore, the stellar response density max¬ 
ima represent the regions of the arms where stars 
would crowd for long timescales, this is, where the 
existence of stable, long-lasting spiral arms would 
be more likely. On the other hand, if the stellar 
density response forms a spiral arm with a differ¬ 
ent pitch than the imposed angle, then the im¬ 
posed spiral arms triggered on the disk (by a bar, 
an interaction, etc.), would rather be structures of 
transien t nature since those are no t orbitally sup¬ 
ported (jPerez-Villegas et al. I 12013 I) . Likewise, for 
the case of the spiral arms in the Milky Way, the 
values frequently seen in literature for the pitch 
angle, range from ~ 11° to 19°. These values and 
the knowledge of the galactic type could provide 
some information about their nature, i.e. whether 
they are long-lasting or transient structures. Fol- 
lo wing the pitch angle restric tions found 
by IPerez- Villegas et alJ (|2Q13l L the larger 
pitch angle values reported in literature for 
the Milky Way galaxy would imply that its 
spiral arms are a transient feature. The for¬ 
mation of four spiral arms in the gas response is, 
in this scenario, another piece of evidence of a 
transient nature of the spiral arms in the Milky 
Way galaxy, as we claim it is the secondary pair 
of arms in Figure [9l that coincides more closely 
with the stellar density re sponse from p eriodi c 
orbits, as expected from iGomez et al.l (|2013l b 


Summarizing, the first (imposed) pair of massive 
spiral arms formed in the disk, with a larger pitch 
angle, triggers a second pair of arms (traced ap¬ 
proximately by the periodic orbits), wit h smaller 
pitch angles (jPerez-Villegas et~al] l2013h . In this 
outline, the gas responds forming a second pair of 
arms aligned with the locus of the orbital crowd¬ 
ing. This lighter structures would likely be pref¬ 
erentially formed by young stars and gas than by 
evolved stars because of their transient nature, 
i.e. similar to what we call “branches”. Finally, in 
this framework, the presence of clear and strong 
branches in spiral galaxies, with smaller pitch an¬ 
gles than the corresponding couple of massive spi¬ 
ral arms on a galaxy, would be a signature of 
the transient nature of the spiral arms in a given 
galaxy. On the other hand, spiral galaxies without 
evidence of branches could indicate the presence 
of a longer lasting spiral arm structure. 


4 DISCUSSION AND CONCLUSIONS 

With the use of a detailed three dimensional, 
density-distribution based potential for the spi¬ 
ral arms, combined with MHD simulations on 
a Milky Way-like galactic disk, we have studied 
the stellar orbital and gaseous response to the 
galactic potential. As a first experiment, we con¬ 
structed a simple cosine potential (as the ones 
commonly employed in literature) that repro¬ 
duced approximately what the density based po¬ 
tential PERLAS exerts on the stellar and gaseous 
dynamics. 

The first set of simulations compare the gas 
response when the disk is perturbed by both 
spiral arm potential models. We also varied the 
structural parameters of the spiral arms in the 
PERLAS model, such as the pitch angle (12° to 
19°) and the mass of spiral arms (1.75% to 5% of 
the stellar disk mass) in order to understand how 
these parameters affect the gaseous disk dynam¬ 
ics. Additionally, we constructed stellar periodic 
orbits and calculate the stellar response density 
maxima. 

With these exercises we found that only in 
the case of the PERLAS model the gas and stel¬ 
lar density response (based on the existence of 
periodic orbits) is a consistent four-armed spiral 
structure: a couple of strong gaseous arms located 
at the position of the imposed stellar arms, and a 
second pair of weaker gaseous arms located at the 
position of the stellar orbit crowding. Since the 
potential and the stellar density response 
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Figure 9. Stellar and gaseous response to different parameters of the spiral arm model (PERLAS). Panels on the 
left column show the stable periodic stellar orbits {black solid lines)^ the response density maxima {filled squares), 
and the imposed spiral arm locus {open squares), flanked by two dotted lines that represent the spiral arm width. 
Colored circles represent the position of the inner Lindblad {green), 4/1 {blue) and corotation {yellow) resonances. 
Panels on the right column show the gas density distribution after 1.2 Gyr, together with the inner Lindblad, 4/1 
and corotation resonances {dashed eircles). The mass of spiral arms and pitch angle are: Mgp = O.OITSM^:) with 
i = 12° {first row) and i = 19° {second row), Msp = 0.05Mjj with i = 12° {third row) and i = 19° {fourth row). 
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do not coincide, the spiral arms are prone 
to destruction. So, the presence of the sec¬ 
ond (gaseous) pair is interpreted as a sign 
of this transiency. 

Our main conclusions can be summarized as 
follows: 

• We performed a study of the gas response 
to two galactic potentials: the density distribu¬ 
tion based model PERLAS and the widely em¬ 
ployed in literature cosine potential. We verified 
that the gas response (a two armed structure) 
is the same for both models close to the lin¬ 
ear regime only, i.e. for low-mass spiral arms 
and pitch angles smaller than ~ 10°. In the 
general case however (i.e. large pitch angle or 
arm mass), even when the cosine potential has 
been fitted as close as possible to the PERLAS 
model, they produce quite distinct outcomes on 
the gas response. In the case of the cosine po¬ 
tential, the gas responds invariably forming two 
spiral arms while, with the PERLAS model, the 
gas responds with four gaseous spiral arms (in the 
general strong arm case). 

• We increased the strength of the spiral arms 
represented by the cosine potential up to a point 
the arms were equivalent and even beyond the 
mass of a strong bar as an experiment to try to 
reproduce the gas response provided by the PER¬ 
LAS model. However, the answer was always a 
bisymmetric structure. We conclude that the spi¬ 
ral arm strength is not responsible for the four- 
arm gas response, but rather it is the product of 
the forcing generated by the whole density distri¬ 
bution better represented by the PERLAS model 
that, in turn, forces the periodic orbit response 
to shift its crowding regions inside the imposed 
locus of the massive spiral arms. 

• Using the PERLAS model, we changed the 
structural parameters of the spiral arms accord¬ 
ing to the observational and theoretical uncer¬ 
tainties in the determination of the Milky Way’s 
spiral arms in literature. Within these values, 
for the general case, the stellar response den¬ 
sity maxima systematically forms spiral arms 
with a smaller pitch angle than the imposed 
spiral, meaning that the spiral arms might be 
a transient feature i n the Milky Way Galaxy 
(|Perez- Villegas et~al] l20I3l L The presence of a 
second pair of lighter spiral arms (“branches”), 
with smaller pitch angles induced by the first 
more massive stellar spiral arms, might be evi¬ 
dence of the lack of support to the stronger arms 
on a galactic disk, and therefore evidence of the 


transient nature of spiral arms in a given galaxy. 
Applying this scenario to the Milky Way, for the 
stronger spiral arm values reported in literature 
(i.e. pitch angles larger than ~ 10°, and masses 
larger than ~ 2% of the disk), the spiral arms in 
the Milky Way would be of transient nature. 

• Although in this work we applied the mod¬ 
els to Milky Way-like galactic discs, it is worth 
noticing that the results are general. This means 
that the presence of branches with smaller pitch 
angles than the main arms might be the “smok¬ 
ing gun” that proves transiency of spiral arms in 
any given galaxy. 


ACKNOWLEDGMENTS 

We thank Edmundo Moreno for enlightening dis¬ 
cussions that helped to improve this work. APV 
acknowledges the support of the postdoctoral Eel- 
lowship of DGAPA-UNAM, Mexico. This work 
has received financial support from DGAPA PA- 
PIIT grants INI 11313 and INII4II4. 


REFERENCES 

Allen, G. V Santillan, A. 1991, Rev. Mexicana 
Astron. Astrofis., 22, 256 
Amaral, L. H., & Lepine, J. R. D. 1997, MNRAS, 
286, 885 

Antoja, T., Eigueras, E., Torra, J., Valenzuela, 
O., & Pichardo, B. 2010, Lecture Notes and 
Essays in Astrophysics, 4, 13 
Antoja, T., Valenzuela, O., Pichardo, B., 
Moreno, E., Eigueras, E. & Eernandez D. 2009, 
ApJ, 700, L78 

Antoja, T., Romero-Gomez, M., Eigueras, E., 
Valenzuela, O., Pichardo, B. & Moreno, E. 
2011, MNRAS, 418, 1423 
Artymowicz, P., V Lubow, S. H. 1992, ApJ, 389, 
129 

Balbus, S. A. 1988, ApJ, 324, 60 
Barbanis, B., V Woltjer, L. 1967, ApJ, 150, 461 
Bate, M. R., 1997, MNRAS, 285, 16 
Benjamin, R. A., Ghurchwell, E., Babler, B. L., 
et al. 2005, ApJ, 630, LI49 
Buta, R. V Block, D. L., 2001, ApJ, 550, 243 
Ghakrabarti, S., Laughlin. G., V Shu, Erank H., 
2003, ApJ, 596, 220 

Gombes, E., V Sanders, R. H. 1981, AVA, 96, 
164 

Gontopoulos, G. 1967, Bull. Astron. (Ser. 3), 2, 
223 


© RAS, MNRAS 000, [IHTT] 






16 Perez-Villegas, Gomez & Pichardo 


Contopoulos, G. & Grosbpl, P., 1986, A&A, 155, 
11 

Gontopoulos, G. & Grosbpl, 1988, A&A, 197, 83 
Gontopoulos, G., & Harsoula, M. 2012, Gelestial 
Mechanics and Dynamical Astronomy, 113, 81 
Gorder, S., Sheth, K., Scoville, N. Z., et ah 2008, 
ApJ, 689, 148 

Dobbs, G. L., k Bonnell, L A. 2006, MNRAS, 
367, 873 

Dobbs, G. L., Burkert, A., k Pringle, J. E. 2011, 
MNRAS, 417, 1318 

Donner, K. J., k Thomasson, M. 1994, A&A, 
290, 785 

D’Onghia, E., Vogelsberger, M., k Hernquist, L. 

2013, ApJ, 766, 34 
Drimmel, R. 2000, AkA, 358, L13 
Drimmel, R., k Spergel, D. N. 2001, ApJ, 556, 
181 

Dwek, E., Arendt, R. G., Hauser, M. G., et ah 
1995, ApJ, 445, 716 

Eden, D. J., Moore, T. J. T., Morgan, L. K., 
Thompson, M. A., k Urquhart, J. S. 2013, MN¬ 
RAS, 431, 1587 

Elmegreen, D. M. 1980, ApJ, 242, 528 
Englmaier, P., k Gerhard, O. 1999, MNRAS, 
304, 512 

Eeitzinger, J.V., k Schwerdtfeger, H., 1982, 
AkA, 116, 117 

Eoyle, K., Rix, H.-W., Walter, E., k Leroy, A. K. 
2010, ApJ, 725, 534 

Eoyle, K., Rix, H.-W., Dobbs, G. L., Leroy, 
A. K., k Walter, E. 2011, ApJ, 735, 101 
Eujimoto, M. 1968, Non-stable Phenomena in 
Galaxies, 29, 453 
Eux, R. 1999, A&A, 345, 787 
Gomez, G. G., Pichardo, B., k Martos, M. A. 
2013, MNRAS, 430, 3010 
Grosbpl, P. J. k Patsis, P. A. 1998, AkA, 336, 
840 

Junqueira, T. G., Lepine, J. R. D., Braga, 
G. A. S., k Barros, D. A. 2013, A&A, 550, 
AA91 

Kawata, D., Hunt, J. A. S., Grand, R. J. J., 
Pasetto, S., k Gropper, M. 2014, MNRAS, 443, 
2757 

Kim, W.-T., Kim, Y., k Kim, J.-G. 2014, ApJ, 
789, 68 

Kim, Y., k Kim, W.-T. 2014, MNRAS, 440, 208 
Kim, W.-T. k Ostriker, E. G., 2002, ApJ, 570, 
132 

Kim, W.-T. k Ostriker, E. G., 2006, ApJ, 646, 
213 

Laurikainen, Eija k Salo, Heikki, 2002, MNRAS, 


337, 1118 

Lee, Wing-Kit k Shu, Erank H., 2012, ApJ, 756, 
45 

Lin, G. G. k Shu, E. H., 1964, ApJ, 140, 646 
Lynds, B. T. 1970, The Spiral Structure of our 
Galaxy, Proc. lAU Symp. 38, 26 
Martinez-Garcia, E. E., k Gonzalez-Lopezlira, 
R. A. 2013, ApJ, 765, 105 
Martos, M., Hernandez, X., Yanez, M., Moreno, 
E., k Pichardo, B. 2004, MNRAS, 350, L47 
Matsumoto, T., Hayakawa, S., Koizumi, H., et 
al. 1982, The Galactic Genter, 83, 48 
Moore, T. J. T., Urquhart, J. S., Morgan, L. K., 
k Thompson, M. A. 2012, MNRAS, 426, 701 
Patsis, P. A., Hiotelis, N., Gontopoulos, G., k 
Grosbol, P. 1994, A&A, 286, 46 
Patsis, P. A., Grosbol, P., k Hiotelis, N. 1997, 
AkA, 323, 762 

Perez-Villegas, A., Pichardo, B., k Moreno, E. 
2013, ApJ, 772, 91 

Perez-Villegas, A., Pichardo, B., Moreno, E., Pe- 
imbert. A., k Velazquez, H. M. 2012, ApJ, 745, 
L14 

Pichardo, B., Martos, M., Moreno, E. k Espre- 
sate, J., 2003, ApJ, 582, 230 
Roberts, W. W. 1969, ApJ, 158, 123 
Roberts, W.W., Jr., Huntley, J. M., k van Al- 
bada, G. D., 1979, ApJ, 233, 67 
Roca-Eabrega, S., Valenzuela, O., Eigueras, E., 
et al. 2013, MNRAS,432, 2878 
Roca-Eabrega, S., Antoja, T., Eigueras, E., et 
al. 2014, MNRAS, 440, 1950 
Roskar, R., Debattista, V. P., Quinn, T. R., k 
Wadsley, J. 2012, MNRAS, 426, 2089 
Sanders, R. H., k Tubbs, A. 1980, AJ, 235, 803 
Scarano, S., V Lepine, J. R. D. 2013, MNRAS, 
428, 625 

Scoville, N. Z., Polletta, M., Ewald, S., et al. 
2001, AJ, 122, 3017 

Sellwood, J. A. 2011, MNRAS, 410, 1637 
Sellwood, J. A., k Garlberg, R. G. 2014, ApJ, 
785, 137 

Seo, W.-Y., k Kim, W.-T. 2014, ApJ, 792, 47 
Shetty, R., k Ostriker, E. G. 2006, ApJ, 647, 997 
Shu, E. H., Milione, V., k Roberts, W. W., Jr. 
1973, ApJ, 183, 819 

Stone J. M., Norman M. L., 1992a, ApJS, 80, 
753 

Stone J. M., Norman M. L., 1992b, ApJS, 80, 
791 

Sofue, Y. 1976, AVA, 48, 1 
Tsoutsis, P., Efthymiopoulos, G., k Voglis, N. 
2008, MNRAS, 387, 1264 


© RAS. MNRAS QOQ.mfTTl 


The Galactic Branches as Evidence for Transient Spiral Arms 


Vallee, J. P. 2002, ApJ, 566, 261 
Vallee, J. P. 2013, International Journal of As¬ 
tronomy and Astrophysics, 3, 20 
Voglis, N., Stavropoulos, L & Kalapotharakos, 
C. 2006, MNRAS, 372, 901 
Wada, K., Baba, J., & Saitoh, T. R. 2011, ApJ, 
735, 1 

Wada, K., & Koda, J. 2004, MNRAS, 349, 270 
Weiner, B. J., & Sellwood, J. A., 1999, ApJ, 524, 
112 

Weaver, H. F. 1970, Interstellar Gas Dynamics, 
lAU Symp. 39, 22 

Yano, T., Kan-Ya, Y., & Gouda, N. 2003, PASJ, 
55, 409 

Yanez, Miguel A., Norman, Michel L, Martos, 
MArco A., & Hayes, John G., 2008, ApJ, 672, 
207 

Zhang, X. 1998, ApJ, 499, 93 
Zhao, H. 1996, MNRAS, 283, 149 


© RAS, MNRAS 000, [IHTT] 


